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Abstract 

In this paper we report on two related topics; Kinetic Monte Carlo simulations of the steady state growth 
of rod eutectics from the melt, and a study of the surface roughness of binary alloys. 

We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair 
exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the 
pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid 
and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction 
is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, 
arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically 
have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt 
model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or non- 
faceted phases depending on the entropy of fusion. 

There have been many studies of the surface roughening transition in single component systems, but none 
for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy 
determines whether the eutectic morphology will be regular or irregular. We have conducted a study of 
surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found 
to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase 
diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness. 

Introduction 

The goal of this project is to gain insight into the kinetics of growth in faceted / non-faceted rod eutectic 
(irregular eutectic) alloys. These eutectic alloys often form irregular microstructures. Eor example, 
controlling the microstructure uniformity in irregular eutectics such as Ge-TiGCj and Si-TaSi, is critical 
in controlling their structural and electronic properties. Regular eutectics, where both phases are non- 
faceted, closely follow the Jackson-Hunt model, whereas irregular eutectics usually grow with a larger 
inter-phase spacing than predicted. This behavior is highlighted with the discrepancies reported between 
ground based and micro-gravity measurements.^ '* 

Irregular eutectics have previously been modeled using adaptations of the Jackson-Hunt model 
however these analytical models do not address the nucleation of new atomic layers on the faceted phase. 
We have conducted kinetic Monte Carlo simulations (kMC) of eutectic growth based on the spin-one 
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Ising model, which provides a definitive description of the solid-liquid interface above and below the 
surface roughening transition^ This model will be used to gain insight into the growth kinetics of systems 
comprised of both faceted and non-faceted solid phases. 

The second part of this report examines the interface roughness of a binary alloy. Since regular and 
irregular eutectics are classified by the surface roughness of the solid phases, understanding the surface 
roughening behavior is important. The location of the surface roughening transition for the two phases 
of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. Below the 
roughening transition, there is nucleation barrier to the formation of new layers exists, and the growth rate 
depends on the kink site density.® '® As a result, the undercooling required for growth on smooth interfaces 
below the transition point is much greater than for interfaces above the surface roughening transition, so 
the two behave quite differently. 


Kinetic Monte Carlo 

Simulations for this study were conducted using a spin-one Ising model with the liquid assumed to be an 
ideal solution in contact with a regular solution solid, where both phases are comprised of A and B atoms. 
Simulations are implemented on an initially fiat solid-liquid interface. Solid atoms are added to the solid 


from the liquid at a rate of; 
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Here AS. is the entropy difference between the liquid and the solid phases of the A and B atoms, is the 
pre-exponential factor with units of 1/time, and k^ is the Boltzmann’s constant. Correspondingly, atoms 
leave the solid and join the liquid at a rate of: 
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Here O.. represents A- A, B-B, and A-B bonds in the solid and T is the temperature of the system. Both 
solid-liquid bonds and liquid-liquid bonds energies are set to zero. Simulations reported for this study are 
all done isothermally. The number of both A and B atoms was conserved in the simulation cell. Diffusion 
was conducted only in the liquid and was implemented by pair exchange of B atoms at a rate proportional 
to add-atom arrival. 

v^=i3v""=— (3) 


where 0 is the proportionality constant. This rate may also be written in terms of the diffusion coefficient 
D and the cube root of the atomic volume a. 
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Figure 1 . Simulation cell geometry. 

Rods of solid B atoms are arranged in the simulation cell with semi-periodic boundary conditions in a 
closed packed formation. This is done with a simulation cell of R by V3 R where R is the rod spacing. 

The rod radius r is dependent on the eutectic composition. 

These kMC simulations were implemented in a similar fashion to previous work," events were chosen one 
at a time from a list of all possible events. The probability of choosing an event was inversely proportional 
to the frequency of that event. After the event is implemented the list is updated. 

Rod Eutectic Simulations 

The kMC simulations of a rod eutectic system started with a simulation cell where a liquid of randomly 
distributed A and B atoms was in contact with a solid consisting of columns of pure B embedded in an A 
matrix. The rods were arranged in a close packed formation in the simulation cell with periodic boundary 
conditions normal to the growth direction, as displayed in figure 1 . The rod positions and radii are not 
confined in the simulation cell, however the simulation cell does place an upper limit on the rod spacing. 
Thus, the size of the simulation cell is an important parameter in these simulations. 
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Figure 2. Sample phase diagram. 

In the binary eutectic systems studied, the a and P phases contained very little of the second component. 
The solidus lines in the phase diagram above are indistinguishable from the axis. This phase diagram 

was derived from the following bond energies: 

/ k3= 2820 “K, AS^ = 2.5 k^, / k^ = 3000 “K, AS^ = 5.0 k^, / k^ = 500 “K 

Phase Diagram 

Equilibrium phase diagrams for our eutectic systems were calculated for an ideal liquid in contact with 
a regular solution solid. In order to examine rod eutectic systems we used phase diagrams where the 
eutectic composition of the minor component was at 10%. The a and P phases in the eutectic systems 
used in this study contained very little of the second component. 

Jackson-Hunt model 

Starting values for the rod spacing R, which depends on the diffusion coefficient D and the simulation 
temperature, were derived from the Jackson-Hunt modeled. The temperature of the simulation, T = T^ 
- AT, was chosen for reasonable rod spacing: 

Where T^ and are the eutectic temperature and composition and m^ and mp are the slopes of the liquidus 
lines of the a and P phases at the eutectic point. Rod spacing R of ~90 unit cells were chosen based on 
computational restraints. 

The rod spacing is related to the diffusion coefficient (D) by: 
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Figure 3. Typical simulation results. 

(a) and (b) are 3D perspective views of the solid atoms in the simulation cell, (c) is a diagonal slice 
though the simulation cell of figure (b) where A and B solid atoms and B liquid atoms are displayed. 
The simulations cell has dimensions of 50 by 87 by 230 unit-cells. Figures were rendered using 

Raster3D.'^ 


®AA/k3=2820“K,AS, = 2.5k3 
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Results 

Simulations of stable growing rod eutectics were based on a simple cublic lattice oriented in the {100} 
growth direction. Stable rod growth was achieved when the liquid was 3% richer in the minor component 
as compared to the equilibrium phase diagram. The rods were also observed to modulate in thickness when 
the liquid was off this composition. Both of these behaviors are experimentally observed characteristics. 
An example of one of these simulations is displayed in figure 3. 

Irregular Eutectics 

Interest in examining the roughening behavior of the binary alloy was motivated by the simulation above. 
With a simple cubic simulation cell oriented in the {111} direction and both A and B atoms having the 
same AS’s, we observed the surface roughness of the two phases was clearly different. The minor phase 
clearly showed faceting while the interface of the major phase was rough. We also observed that it was more 
difficult to get stable growth of faceted rods if the facets were in the plane of the growth direction. At this 
point we changed our focus from SC to FCC lattices since the SC {111} simulations cell does not tessellate 
in the same way as the SC {100} simulation cell. Since FCC {111} planes have more in the plane nearest 
neighbors than { 100 } direction systems maybe rough { 100 } faces while facets in the { 111 } direction. 
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(a) (b) 

Figure 4. Irregular rod eutectic. 

Simulation conducted on a simple cubic lattice oriented in the {111} growth direction. Both A and B 
atoms gave the same Delta S however the surface roughness characteristics of both the alpha and beta 
phase are clearly different. Faceting is along of the {100} faces. The 3D perspective view of figure (a) 
is distorted into a rectangular simulation cell. Figure (c) is a diagonal slice thought the simulation cell. 

2000 "K, AS, = 5.0 k^, / k = 1650 "K, AS^ = 5.0 k„, / k„ = 300 “K 




1028.2 “K, AT^ 


9.2 “K 


Surface Roughness 

While the surface roughening transition of single component systems has been studied, the surface 
roughening transition for binary systems has not. We have conducted a study of the inter-facial roughness 
in alloys using the spin-one Ising Model with diffusion using kMC to understand the roughening behavior 
of the individual phases in the rod eutectic system. We have found that the surface roughness depends on 
the bond energies and melting temperature as given by the liquidus temperature on the equilibrium phase 
diagram. 


The surface roughness can be defined in terms of the density of lateral bonds missing at the surface. In this 
study, the definition of surface roughness which was used is the same as used in BCF'®: 




( 6 ) 


Here Ng^(t) is the number of solid-liquid bonds on the surface and Ng^(O) is the number of solid-liquid 
bonds on an atomically smooth surface. This does not define the location of the surface roughing transition, 
since the surface roughness increases smoothly through the transition. However, the location of the surface 
roughening transition can be determined from height-height correlations in the configurations of the 
interface. 


For melt growth of single component materials, the location of the surface roughening transition is given 
by the Jackson a-factor. 
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Here L is the latent heat of fusion, T^^is the melting temperature of the system, r] is the number of bonds 
in the plane of the interfaee, and Z is the coordination number of the lattice. 

Surface Roughness Simulations 

The relationship between surface roughness and the a-factor was examined on a FCC lattice oriented in 
the {100} direction. Simulation cells typically had dimensions of 200 by 200 by 100 unit cells. Solid- 
liquid interfaces were created flat with pure solid of A atoms and a mixture of A and B atoms in the liquid. 
The evolution of structure of the interface then proceeded at the liquidus temperature of the binary alloy. 

Surface Roughness Results 

Three separate phase diagrams using different entropies of fusion for the pure a phase were used in the 
simulations reported here. System 1 was run with a AS of 7.0 k^, system 2 with a AS of 6.0 and system 

3 with a AS of 5.0 k^. Pure systems were also simulated over a similar range of a-factor values. Surface 
roughness as a function of a-factor values is reported in flgure 5. Rather than using the melting point of 
the pure material to calculate the a-factor for the alloys, the liquidus temperature of the binary systems 
was used. Error bars on the data are the standard deviation of the fluctuations of the roughness value 
around the reported mean. Different a-factors values apply for each point along the liquidus lines of each 
phase diagram. From flgure 5 it is evident that both pure and binary systems follow the same roughness 
relationship with a smooth monotonically decreasing roughness with increasing a-factor. 
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Figure 5. a-factor vs. roughness parameter 
Both binary and pure systems follow the same roughness relationship 
System 1: AS^^^^ = 7.0 k^, = 1500$ “K, = 2000$ °K, = 500$ “K 

System 2: AS^^J= 6.0 k^ “K, = 1470 $ °K, = 2000$ “K, = 500$ “K 

System 3: AS^^^^ = 5.0 k^, = 1310$ “K, = 1800$ “K, = 500$ “K 
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In the case of this simple eutectic system the surface roughness of both phases increases as the system 
moves to either side of the eutectic. In a rod eutectic system this is more dramatic for the minor phase 
because of the larger difference between the melting temperature of the pure material and the alloy. 

Conclusions 

We have developed and demonstrated a kMC spin one Ising model capable of simulating the growth of 
stable rod eutectics. The simulations use the Jackson-Hunt model as a guide to provide appropriate starting 
conditions. These simulations exhibit behaviors similar to experimentally observed eutectic structures. 
While irregular rod eutectic growth has been more challenging, simulations where the minor phase is 
faceted have been more successful when the faceted interface has been out of the plane of the growth 
direction. Since the surface roughness of the two solid phases plays an important role in the development 
of the microstructure, an investigation of the surface roughness of binary alloys was conducted. It was 
found that the surface roughness of alloys follows the same behavior as pure materials provide that the 
liquidus temperature is used to calculate the a-factor. The entropy difference between pure solid and pure 
liquid of the same material is not the same as the entropy difference between the solid and the liquid alloy, 
and this difference is expressed in terms of the lowered melting point. The surface roughening of a phase 
in a eutectic alloy does not depend on the entropy of fusion of the pure components: the interfaces of the 
alloy are smoother than the interfaces of the pure materials. 
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